Fjords in viscous fingering: Selection of width and opening angle 
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Our experiments on viscous fingering of air into oil contained between closely spaced plates reveal 
two selection rules for the fjords of oil that separate fingers of air. (Fjords are the building blocks of 
solutions of the zero-surface-tension Laplacian growth equation.) Experiments in rectangular and 
circular geometries yield fjords with base widths X c /2, where A c is the most unstable wavelength 
from a linear stability analysis. Further, fjords open at an angle of 8.0° ± 1.0°. These selection rules 
hold for a wide range of pumping rates and fjord lengths, widths, and directions. 

PACS numbers: 47.54.-r, 47.20.Ma, 68.35.Ja 



Similar growth patterns have been found for viscous 
fingering, bacterial growth, flame propagation, dielectric 
breakdown, electrodeposition, solidification, metal corro- 
sion, and diffusion- limited aggregation (DLA) 0. Analy- 
ses of these interfacial patterns is daunting because of the 
broad range of length scales, as illustrated by our exper- 
iments on fingering of air in thin oil layers fFig.lTjI. The 
interface becomes increasingly ramified with increasing 
forcing strength, which is characterized by the capillary 
number 



Ca: 



(1) 



where /i is the viscosity of the oil, V is the local inter- 
facial velocity, and a is the surface tension of the air-oil 
interface. 

Most theoretical work on viscous fingering has con- 
centrated on fingers rather than fjords, but fjords have 
been considered as building blocks of a theory of Lapla- 
cian growth in an analysis of a two-dimensional invis- 
cid fluid penetrating a viscous fluid (with zero interfacial 
tension) Q. The exact solutions for the interface are 
free of finite-time singularities and are linear combina- 
tions of logarithmic terms, each term representing a sin- 
gle straight fjord with parallel walls. While these analytic 
solutions have helped in the understanding of the growth 
of complex growth patterns, experiments indicate fjords 
can have curved trajectories and nearly always have a 
nonzero opening angle (Fig.QJ 0- 

Most experiments have also concentrated on fingers, 
but a few experiments have examined fjords For the 
circular sector ("wedge") geometry, Thome et al. mea- 
sured the angular gap between a wall and a divergent 
finger (at a point 12 cm from the vertex) to be approx- 
imately 5° for divergent sector angles between 20° and 
50° p|. Their sector decomposition method considered 
finger growth in terms of wedges whose "virtual walls" 
were fjord centerlines; this approach suggests that a fjord 
opening angle should be about 10°, i.e. twice the angjte 
measured by Thome et al. Q . Theoretical analysis 0, 



yielded an estimate corresponding to a fjord opening an- 
gle of 11.7° @. Other work has suggested that the width 
of a fjord relative to the width of the splitting finger has 
values between | and | [S 0, ■ 

Despite much effort, selection laws are lacking for ram- 
ified fingers in either rectangular [Figs.^a) and Qlb)] or 
circular geometries [Figs, Hfc) and [T^d)] . We find that 
fjords in both geometries have a well-defined base width 
that is determined by the characteristics of the interface 
at the fjord's birth (when a finger splits). Further, we 
find that fjords open at an angle that is close to 8°. The 
selection rules for the fjord width and opening angle hold 
for a wide variation in the ramification of the patterns. 

Experimental apparatus — Interfacial patterns were 
grown in both rectangular and circular Hele-Shaw cells 
(closely spaced glass plates) filled with oil. Air was forced 
into the oil, creating an unstable interface (Fig. [TJ. in 
the rectangular geometry, two pieces of glass (each 1.91 
cm thick), separated by a gap of 0.0508 cm, bounded a 
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FIG. 1: Viscous fingering patterns of air (black) penetrating 
into oil (white) in rectangular and circular Hele-Shaw cells for 
different forcing levels. A fjord is the region of oil between two 
adjacent fingers. Rectangular cell: Ca = (a) 0.0056 and (b) 
0.046; the rectangles are 25 cm wide and 190 cm long, while 
the entire cell is 254 cm long. Circular cell, Ca =(c) 0.000 88 
and (d) 0.27; the radial patterns are approximately 20 cm in 
diameter, grown in a cell of diameter 28.8 cm. The Ca values 
were computed using the average velocity of the farthest tip. 




FIG. 2: (a) Time development of a fjord grown in the rectangular cell at Ca = 0.0076 (cf. Fig. 1 caption). Adjacent interfaces 
are separated by 1 s. (b) Curvature k and (c) velocity V of the point on the finger (k > 0) that evolves into a fjord (k < 0) 
and stops moving; t = is the time at which « = 0. The velocity Vo used in computing the instability wavelength A c [Eq. 
corresponds to k = 0; here Vo = 0.25 ± 0.01 cm/s. The arrow points to the minimum of curvature in time, which defines the 
fjord centerline arc length £ — 0. (d) Fjord width Ifasa function of the fjord centerline arc length I. The base width of the 
fjord, Wo, is given by the fjord width extrapolated to I — 0; here Wo = 1.08 ± 0.06 cm. (e) The selected opening angle of the 
fjord, given by the plateau in tp(£), is ?/>o = 7.1° ± 0.5°. 



cell 25.4 cm wide by 254 cm long 9j . The maximum vari- 
ation in gap thickness at the center of the cell (due to the 
imposed pressure gradient) was less than 3% even at the 
highest pump rate . The cell was filled with silicone oil 
(viscosity /i = 50.6 mPa s and surface tension a = 20.5 
mN/m at 24 °C); the oil wets the glass completely. A 
uniform flow rate was achieved using a syringe pump to 
withdraw oil from a reservoir at one end of the channel; 
an air reservoir at atmospheric pressure was attached to 
the other end. The channel was illuminated from below. 
Each full interface, as in Figs, ^a) and ^b), was con- 
structed from 11 images (each 1300 x 1030 pixels, 0.21 
mm/pixel), which were obtained using a camera and ro- 
tating mirror. Alternatively, the camera was focused on 
a fixed 22-cm-long section of the channel to obtain a time 
resolution of 12 frames/s [Fig. Ela)]. 

In the radial geometry, air was forced into an oil-filled 
gap through a hole in the center of the bottom glass 
plate. Each optically polished glass plate had a diameter 
of 28.8 cm and a thickness of 6.0 cm. The two plates were 
separated by a gap thickness of 0.0127 cm (uniform to 

i%) naim. The gap and an annular reservoir were filled 

with silicone oil (viscosity /i = 345 mPa s and surface 
tension a — 20.9 mN/m at 24 °C). Interfacial patterns 
were grown either by maintaining a constant pressure 
difference between the oil reservoir and air or by using 
a syringe pump to remove oil from the buffer. Images 
of resolution 0.32 mm/pixel were acquired at up to 12 
frames/s. 

Fjord characteristics — From the images we deter- 
mined, as a function of time, the interface's velocity V 
and curvature k, the fjord width W , and the fjord open- 
ing angle ip (Fig. [21 . The most unstable (critical) wave- 
length of a curved interface between a less viscous fluid 



forced into a more viscous fluid is given by 0, Il3] | 



^Ca/6 2 + K 2 /12 



(2) 



where b is the cell gap thickness and k is the interface 
curvature (not including the curvature in the small gap 
dimension) at t he p oint where the finger becomes unsta- 
ble 0, flEl ll6L Il7l |. For all interfaces that we examined, 
the curvature term in Eq. was about two orders of 
magnitude smaller than the term involving the capillary 
number Ca; hence we made the flat interface (n = 0) 
approximation in computing A c . 

When a finger became unstable, the finger flattened 
and a fjord was born. A normal projection algorithm 
was used to track the fjord trajectory back in time to 
earlier interfaces [solid, nearly straight line on the left 
in Fig.|3[a)]. (The fjord location for an earlier interface 
was identified by projecting a ray normally from the later 
interface until it intersected the earlier one.) The time 
evolution of the local curvature k then revealed the birth 
of the fjord, where k = [Fig. Gib)]. The velocity Vq at 
this point [Fig. EJc)] was used in Eq. J5J to calculate \ r . 
The uncertainty in A c is 9%; film wetting corrections [Ij, 
fTEj l and the velocity measurements are the main sources 
of uncertainty. 

Once a fjord developed, we determined its width W as 
a function of the arc length I of its centerline, as illus- 
trated in Fig. GId). The fjord centerline was found by 
beginning with a point on the side of the fjord; then the 
closest point on the opposite side was found such that 
the interior angles between the two local tangents and 
the line connecting the two points were the same. The 
midpoint of the connecting line gave a point on the fjord 
centerline. Repeating this procedure for each point on 
the side of the fjord constructed the centerline. The fjord 
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FIG. 3: The measured fjord base width W [Fig. Eld)] versus 
the predicted width, A c /2, where A c is the instability wave- 
length, Eq. The inset illustrates how an initial infinites- 
imal perturbation of a front sets the wavelength for a large- 
amplitude fjord. Here, the front is moving to the right and is 
observed in a comoving frame. The lengths of the error bars 
are the rms deviations from the line. 



base point (£ = 0) was taken to correspond to the mini- 
mum of curvature along the fjord's trajectory [Fig.|2Jb)]. 
As surface tension relaxed the interface, a fjord's base 
moved forward. The distance a base moved after t = 
was corrected by adding the same distance to the arc 
length [Fig.Efd)]. 

Curves W(£) obtained for different fjords exhibited the 
same structure: W(£) initially grew quickly with increas- 
ing t due to the rounded nature of the fjord base, but soon 
W(£) increased linearly with I. The fjord base width Wq 
was measured as follows: for each £ the local tangent 
of W(£) was linearly extrapolated back to t = 0, a his- 
togram of extrapolated base widths was made, and the 
histogram's peak gave Wq. For an individual fjord mea- 
sured over a wide range of times, this procedure yielded a 
fjord base width with a typical uncertainty of 6%. Con- 
finement by neighboring fjords caused a slightly increas- 
ing extrapolated fjord base width with increasing values 
of L 

The observed linear increase in W(£) suggested an al- 
most constant opening angle: ip = arctan(^i). The an- 
gle tp(£) had the same qualitative evolution for all fjords: 
a fjord began with an opening angle of 180°, followed by 
a surface-tension-induced dip, followed by a peak and a 
plateau (usually with a slightly negative slope), and fi- 
nally a rapid increase at the fjord's end [Fig. Ete)]. We 
defined the opening angle ipo to be the peak of the his- 
togram of ip values generated from plots of ip(£). The 
angle selected was always about 8°; the uncertainty in 
the opening angle for a single fjord was 0.5°. 

Fjord base width — We predict that the base width of 



a fjord is determined at the onset of the instability of the 
moving interface, as illustrated by the inset of Fig.|3 The 
critical wavelength A c given by a linear stability analy- 
sis of the infinitesimally perturbed flat front, Eq. J2J), 
determines the distance between adjacent fingers when 
they are born, and the base width of the macroscopic 
fjord that develops between the two emergent fingers is 
Wq = iA c - A comparison between measurements and 
the predicted widths requires a value for the velocity V 
in the expression for the capillary number Ca. Lajeunesse 
and Couder suggest that a perturbation on a finger is the 
"precursor" of tip splitting and confirm that the pertur- 
bation's wavelength is "of the order of 1.5 A c 3]." Our 
approach tests this idea at the instability's onset, making 
the assumption that V = Vo is measured at the point on 
the fjord's trajectory corresponding to k = [Figs.Hfb) 
and EIc)]. 

Measurements of the fjord base widths Wq in both the 
circular and rectangular geometries agree well with the 
predicted widths, as Fig. OH illustrates. The fjord widths 
varied by an order of magnitude as the control conditions 
were varied widely (pump rates were 0.17 - 2.00 cm 3 /s 
in the rectangular cell and 0.000 83 - 0.033 cm 3 /s in the 
circular cell, and in the latter geometry, measurements 
were also made for pressure differences of 0.10 and 0.25 
atm). The local forcing strength at the birth of a fjord is 
affected by nearby fjords, walls, and details of the inter- 
face, so a given global pumping rate can yield fjords with 
different properties. The uncertainties in A c and Wq are 
9% and 6%, respectively. The observed scatter in Fig. 3 
is consistent with these estimated errors. 

Fjord opening angle — The opening angle -i/>o mea- 
sured in the rectangular geometry is shown in Fig. Ufa) 
as a function of the appropriate forcing parameter for 
this geometry, the modified capillary number 1/B = 
12Ca(w/b) 2 (where w is the channel width) [U^^], which 
was computed using the average velocity of the farthest 
finger tip. Isolated fjords occur in the regime just beyond 
the onset of tip splitting, which occurs for 1/B ss 4000 

MM 

For 1/B > 25 000, the pattern develops a high density 
of fjords. Despite the wide range of experimental param- 
eters in this regime, the measured angles are surprisingly 
similar. The results for the distribution of opening an- 
gles in the circular cell are indistinguishable from those 
for ramified patterns in the rectangular cell: the proba- 
bility distribution function p(ipo) has the same form with 
essentially the same mean 8° and same standard devia- 
tion 1° [Fig. Efb)]. We suggest that the near constancy 
of the opening angle is fundamental to scale invariance 
and fractal dimension. Our result is consistent both with 
efforts to relate the multifractal dimensions of a fully de- 
veloped fingering structure [2(J with the fjords' opening 
angle [2l| and with a recent observation of an invari- 
ant unscreened angle distribution in a fractal grown in 
a radial geometry. While intriguing, the connection be- 
tween these invariant geometrical characteristics remains 
obscure. 

Conclusions — We have presented two selection rules 
for the geometric form of fjords in rectangular and circu- 
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lar cells. The base width selected by the fjords has been 
found to be well described by Wo = iA c , where A c is 
the local critical wavelength. This selection rule provides 
insight into the most apparent difference in viscous fin- 
gering patterns: highly forced patterns are composed of 
thinner fjords. The rule is based on universal features of 
Laplacian growth and does not invoke any specific prop- 
erties of viscous fingering in a Hele-Shaw cell. Therefore, 
we conjecture that this rule should apply not only to 
viscous fingering, but also to other isotropic Laplacian 
growth problems. 



tpo (deg) 



The second selection rule is obtained from measure- 
ments of the opening angle as a function of arc length 
along the fjord centerline: for a broad range of experi- 
mental parameters and for both circular and rectangu- 
lar geometries, the selected opening angle ipo was always 
near 8°; we have no explanation for this result. The 
two selection rules hold for a wide range of fjord widths, 
lengths, and degrees of bending, and for a wide range of 
angles of fjord directions with respect to the rectangular 
cell axis or with respect to radial lines in the circular cell. 



FIG. 4: (a) The fjord opening angle ipo in the rectangular 
cell increases slightly with forcing parameter (the modified 
capillary number 1/B). (b) Probability distributions p(^o) 
measured for the opening angle of all fjords observed in the 
circular cell (•) and of fjords in the high-fjord-density regime 
(1/B > 25 000) of the rectangular cell (A). The data for 
the rectangular and circular cells yield, respectively, ipo = 
7.9° ± 0.8° and ipo = 8.2° ± 1.1°, where the uncertainties are 
one standard deviation. 
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